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One-dimensional scattering problems are of wide physical interest and are encountered in many diverse applications. In this 
article I establish some very general bounds for reflection and transmission coefficients for one-dimensional potential scattering. 
Equivalently, these results may be phrased as general bounds on the Bogolubov coefficients, or statements about the transfer 
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effectively than the corresponding smoothed potentials. The analysis also serves to collect together and unify what would 
otherwise appear to be quite unrelated results. 
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I. INTRODUCTION 

One-dimensional scattering problems occur in a wide 
variety of physical contexts. In acoustics one might be 
interested in the propagation of sound waves down a long 
pipe, while in electromagnetism one might be interested 
in the physics of wave-guides. In quantum physics the 
canonical examples are barrier penetration and reflec- 
tion, while in classical physics an equivalent problem is 
the analysis of parametric resonances. All of these phys- 
ical problems can be analyzed in the same mathematical 
framework, though for definiteness I shall present the dis- 
cussion in terms of the Schrodinger equation, comment- 
ing on alternative formulations as appropriate. 

For one-dimensional scattering problems there is a 
large catalog of specific potentials for which exact an- 
alytic results are known. There are also well-developed 
numerical techniques for estimating the scattering prop- 
erties. In this paper I wish to take a different tack: I 
shall develop a number of very general and rather simple 
bounds on the reflection and transmission probabilities 
(equivalently, these bounds can be presented in terms 
of the Bogolubov coefficients, or in terms of statements 
about the transfer matrix). These bounds, because they 
are so general, are powerful aids in the qualitative un- 
derstanding of one-dimensional scattering. Furthermore, 
this analysis provides a unifying theme that serves to con- 
nect together seemingly quite disparate results obtained 
in individual special cases. 



II. GENERAL ANALYSIS 



If the potential asymptotes to a constant, 



(2) 



then in each of the two asymptotic regions there are two 
independent solutions to the Schrodinger equation 



±oo; x) 



exp(±ik±oox) 



Vk 



(3) 



Here the ±i distinguishes right-moving modes e +lkx from 
left-moving modes e~ lkx , while the ±oo specifies which 
of the asymptotic regions we are in. Furthermore 



v ±oo 



y/2m(E - V ±OD ) 



(4) 



To even begin to set up a scattering problem the mini- 
mum requirements are that the potential asymptote to 
some constant, and this assumption will be made hence- 
forth. 

The so-called Jost solutions jl(| are exact solutions 
J± (x) of the Schrodinger equation that satisfy 



J+(x 



+oo 



exp(+ik +oc x) 



J_ (x — > — oo) 



exp(-ik^oox) 



(•5) 



(6) 



and 



exp(+i/c_ 00 a;) ex^-ik-mx) 
J+(x — > — oo) — > a -j= h p -j= , (7) 



A. Shabat— Zakharov systems 

Consider the one-dimensional time-independent 
Schrodinger equation [1-15] 



2m dx 2 



ip(x) +V(x) ijj(x) =Eip(x). 



(1) 



„exp(-z/c +oc ,a;) exp(+ifc +oc a;) 
J-[x — > +oo) — > a j== h p 



Vk 



(8) 



Here a and are the (right-moving) Bogolubov coeffi- 
cients, which are related to the (right-moving) reflection 
and transmission amplitudes by 



1 



a ' a 



(9) 



These conventions correspond to an incoming flux of 
right-moving particles (incident from the left) being par- 
tially transmitted and partially scattered. The left- 
moving Bogolubov coefficients are just the complex con- 
jugates of the right-moving coefficients, however it should 
be borne in mind that the phases of j3 and /3* are physi- 
cally meaningless in that they can be arbitrarily changed 
simply by moving the origin of coordinates. The phases 
of a and a* on the other hand do contain real physical 
information. 

In this article I will derive some very general bounds 
on \a\ and \f3\, which also lead to general bounds on the 
reflection and transmission probabilities 



R 



T 



(10) 



The key idea is to re-write the second-order Schrodinger 
equation as a particular type of Shabat-Zakharov Jl7j ] 
system: a particular set of two coupled first-order dif- 
ferential equations for which bounds can be easily estab- 
lished. A similar representation of the Schrodinger equa- 
tion is briefly discussed by Peirls Jlq] and related repre- 
sentations are well-known, often being used used without 
giving an explicit reference (see e.g. |19f| ). However an ex- 
haustive search has not uncovered prior use of the partic- 
ular representation of this article, nor the idea of using 
the representation to place bounds on one-dimensional 
scattering. 

I start by introducing an arbitrary auxiliary function 
p(x) which may be either real or complex, though I do 
demand that <p'(x) ^ 0, and then defining 



ip(x) = a(x) 



i iff 



b(x) 



exp(-ip) 



(11) 



This representation effectively seeks to use quantities re- 
sembling the "phase integral" wavefunctions as a basis 
for the true wavefunction |20[ . This representation is of 
course highly redundant, since one complex number tp(x) 
has been traded for two complex numbers a{x) and b{x) 
plus an essentially arbitrary auxiliary function p{x). In 
order for this representation to be most useful it is best 
to arrange things so that a{x) and b(x) asymptote to con- 
stants at spatial infinity, which we shall soon see implies 
that we should pick the auxiliary function to satisfy 



ip'(x) — * fc± c 



as 



±oo. 



(12) 



To trim down the number of degrees of freedom it is 
useful to impose a "gauge condition" 



d_ 

dx 



d_ 

dx 



'ip / CLX \ V V 

Subject to this gauge condition, 
dip 



= 0. 



(13) 



— = i v ip 1 {a(x) exp(+j(^) — b(x) exp(— ip)} . (14) 
dx 



I now re-write the Schrodinger equation in terms of 
two coupled first-order differential equations for these 
position-dependent Bogolubov coefficients. To do this 
note that 



(Pip _ d 
dx 2 dx 



^i^={ae +i<i = -be- iv }\ (15) 



{llfi T{ae + ^ + be-^} 



VP' 



dx 



;V{ ae +i<P-b e -i<P} 



I (ft 



+i<p _ __ 
dx 



-zip 



tip 



= {ae +lip + be- lv } 



f 

2iip'da c+iv .ip" bc _ ilp 
\/p f dx \fp' 
w' 2 

--^={ae +lv + be- tv } 



up 

2ip> db _ ltp . p" +i 
<p' dx yjipx 



(16) 



(17) 



(18) 



(The last two relations use the "gauge condition".) Now 
insert these formulae into the Schrodinger equation writ- 
ten in the form 



d 2 ip 



dx 2 
to deduce 



— £ = -k( x y ip = - 



2m{E~V{x)) 



(19) 



da 



1 



dx ^W'beM-^) 



+i [k 2 (x) - (p') 2 ] (a + bexp(-2ip)) \, (20) 

db 1 „ . f . 

— =+— ^ a exp(+ 2ip) 

-i [k 2 {x) - {p') 2 } (b + a exp(+2i<p)) j . (21) 

It is easy to verify that this first-order system is compati- 
ble with the "gauge condition" ( |l3| ) , and that by iterating 
the system twice (subject to this gauge condition) one re- 
covers exactly the original Schrodinger equation. These 
equations hold for arbitrary p, real or complex, and when 
written in matrix form, exhibit a deep connection with 
the transfer matrix formalism pl| . 



B. Bounds 

To obtain our bounds on the Bogolubov coefficients 
we start by restricting attention to the case that p(x) 



2 



is a real function of x. (Since ip is an essentially arbi- 
trary auxiliary function this is not a particularly restric- 
tive condition). Under this assumption the probability 
current is 



(22) 



Now at x ~ +00 the wavefunction is purely right-moving 
and normalized to 1, because we are considering one- 
dimensional Jost solutions jl6| . Then for all x we have a 
conserved quantity 



|a| a -|6| a = l. 



(23) 



It is this result that makes it useful to interpret a(x) and 
b(x) as position- dependent Bogolubov coefficients relative 
to the auxiliary function <p(x). Now use the fact that 



d\a\ 



1 



. da da* 



dx 2\a\ V dx dx 



(24) 



and use equation ( |20| ) to obtain 

+i[k 2 - (ip') 2 } [a*bexp(-2itp) - ab* exp(+2i<^)] ) . 

(25) 

That is 

[a*6exp(— 2ip)] ^ . 

(26) 

The right hand side can now be bounded from above, by 
systematically using Ke(A B) < \A\ \B\. This leads to 



d\a\ yW + [fc 2 - (yQ 2 



dx 



2|^| 



\b\. 



(27) 



It is essential that ip be real to have | exp(— 2i(p)\ = 1 
which is the other key ingredient above. Now define the 
non-negative quantity 



%(*),*(*)] = 



\l(<p") 2 + \k 2 {x)-(y>y 
W\ 



and use the conservation law ( p3| ) to write 
d\a\ 



dx 



<$y/\aT^i. 



(28) 



(29) 



Integrate this inequality 



{cosh" 1 |a|} |*' 



< 



dx. 



Taking limits as Xi — > —00 and Xf 



-00 



That is 



cosh 1 \a\ <■ 



\a\ < cosh 



+00 



■d dx. 



<d dx 



Which automatically implies 
< sinh 



dx 



(30) 



(31) 



(32) 



(33) 



Since this result holds for all real choices of the auxiliary 
function tp(x), (subject only to ip' 7^ and p' — > k± oc 
as x — ► ±00), it encodes an enormously wide class of 
bounds on the Bogolubov coefficients. When translated 
to reflection and transmission coefficients the equivalent 
statements are 



and 



T > secli 



R < tanh 



dx 



§dx 



(34) 



(35) 



I shall soon turn this general result into more specific 
theorems. 



C. Transfer matrix representation 



The system of equations (f20[)-(|2l|) can also be written 
in matrix form. It is convenient to define 



(36) 



Then 



d 


a 


1 


ilmjp] 


pexp(— 2ip) 




a 


dx 


b 


~ v 


p* exp(+2ip) 


~ilm[p] 




b 



This has the formal solution 

E(xf,Xi) 



a(x f ) 



a(xi) 
b(xi) 



(37) 
(38) 



in terms of a generalized position-dependent "transfer 
matrix" lEltl 



E(xf,Xi) = 



ilm[p] pexp(— 2ip) 
p* exp(+2iip) —ilm[p] 



dx 



(39) 



3 



where the symbol V denotes "path ordering". In par- 
ticular, if we take Xi — > — oo and xj — * +00 we obtain 
a formal but exact expression for the Bogolubov coeffi- 
cients 



a (3* 
(3 a* 

Pexp 



E(oc, —00) 



1 

v 



ilm[p] pexp(— 2i(p) 
p* exp(+2i<p) — ilm[p] 



dx 



(40) 



The matrix E is not unitary, though it does have deter- 
minant 1. It is in fact an element of the group SU(1, 1). 
Taking 



+1 
-1 



(41) 



then (a z ) 2 = +/, and defining = (E*) T , it is easy to 
see 



E*a z E = a z . 



(42) 



This is the analog of the invariance of the Minkowski met- 
ric for Lorentz transformations in S*0(3, 1). Similarly, if 
we define the "complex structure" J by 



1 

-1 



then J 2 — —I and 



£ f = JEJ. 



(43) 



(44) 



III. SPECIAL CASE 1 

Suppose now that the potential satisfies V+oc — V-oo. 
Also, choose the phase function ip(x) to be tp = koo x. 
We also require k^ ^ 0, that is E > V± OC] . This is the 
special case discussed in a different context by Peirls [jl8) . 
Then the evolution equations simplify tremendously, and 
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\k 2 ~kU _ m\V(x)-V 00 \ 



1k r 



h 2 k„ 



(45) 



Using (hkoo) 2 — 2m(E — Voo), the bounds become 



and 



R < tanh 2 ( - 



ny 2(e-v 00 ) 



\V-Vn\dx), (46) 



[ +0 ° iV-V^dx) . (47) 
J — 00 / 



These bounds are exact non-perturbative results, how- 
ever for high energies it may be convenient to use 



the slightly less restrictive (but analytically much more 
tractable) bounds 



T > 



( f+ c 

TO 

1 _ v 



W-v 



and 



R < 



(f +C 



2Eh A 



\V-Voo\dx 



2Eh z 



(48) 



(49) 



This version of the bounds also holds for all energies, but 
is not very restrictive for low energy. 

The transfer matrices can be analyzed by checking that 
the evolution equations simplify to 



da 
dx 
db 
dx 



-imjV-Vn) 

ft koo 
-im(V - V x ) 

li kao 



{a + 6exp(— 2ikoox)} , (50) 



{a exp(+2ik ocl x) + b} . (51) 



This can be written in matrix form as 

d I" a,l —im(V - Voo) 

dx~ [b\ ~ fi^kZ 

1 exp(— 2ikocx) 

- exp(+2ik oc x) —1 





a 




b 



(52) 



This version of the Shabat-Zakharov system p7| has a 
formal solution in terms of the transfer matrix 



E(xf,Xi) 
V exp 



h 2 k 



Xi 



(V(X) - Voo) 



OO J X f 

1 



-1 



dx 



(53) 



The formal but exact expression for the Bogolubov coef- 
ficients is now 



a (3* 
(3 a* 

Pexp 



= E(po, —00) = 

(V(x) - Voo) 



OO «/ — CXD 



2ik oc x 



1 e 

+2ik oa x _j 



dx 



(54) 



Furthermore, the form of the system (|50|)-(|5l|) suggests 
that it might be useful to define 



a = a exp 
b = b exp 



+ 



fi 2 k, 



00 J — 00 



(V(y) - Voo)dy 
(V(y) ~ Voo)dy 



OO «/ — OO 



(55) 
(56) 
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Then 

da -im(V{x)-V O0 ) ~ 

— = -3- b exp(-2zfc 00 a;), (57) 

ax h fcoo 

db +im(V(x) — I4o) ~ , „ , x /„„x 

— = — d exp(+2ik QC x). (58) 

dx K koo 

This representation simplifies some of the results, for in- 
stance 



a (3* 
j3 a* 

V exp 



h 2 k 



E(oo, — oo) = 

(V(X) - Voo) 

-OO 





oo J — oo 







(59) 



This can be used as the basis of an approximation scheme 
for f3. Suppose that for all x we have \b(x)\ <C 1, so that 
d(x)\ w 1. Then 



dx 



+im(V(x) - V4o) 



exp(+2ifcooa;) 



(60) 



This may be immediately integrated to yield 



h 2 k, 



(V(x) - Voo) exp(+2i/c oc x) dx. (61) 



oo •/ — oo 



This is immediately recognizable as the (first) Born ap- 
proximation. If we instead work in terms of the original 
definition (3 



+im 



h 2 k r 



exp 



+ 



h 2 k, 



+oo 



oo */ — oo 



(V(x) -Voo)dx 
(V(x) — Voo) exp(+2i 

(V(y)-V 00 )dy 



exp 



h 2 k, 



oo J — oo 



dx. 

(62) 



This is one form of the distorted Born wave approxima- 
tion. 

In short, this type of analysis collects together a large 
number of results that otherwise appear quite unrelated. 
By taking further specific cases of these bounds and re- 
lated results it is possible to reproduce many analytically 
known results, such as for delta function potentials, dou- 
ble delta function potentials, square wells, and sech 2 po- 
tentials, as discussed later in this article. (See Section 
VI.) 



IV. SPECIAL CASE 2 

Suppose now we take k(x) = <p'(x). This means that 
we are choosing our auxiliary function so that we use 



the WKB approximation for the true wavefunction as a 
"basis" for calculating the Bogolubov coefficients. This 
choice is perfectly capable of handling the case V+oo 7^ 
V_oo but because of the assumed reality of ip is limited to 
considering scattering over the potential barrier. (This 
is the special case implicit in a different context in |jj"9|] ). 
The evolution equations again simplify tremendously, to 
yield 



da 


1 


dx 


f v 


db 


1 


dx 


f v 



tp" b e^cp(~2iip) 
ip" a exp(+2i<p) 



(63) 
(64) 



This form of the evolution equations can be related to 
the qualitative discussion of scattering over a potential 
barrier presented by Migdal 22 2^|. For this choice of 
auxiliary function 







2\p'\ 

and the bounds become 



2\k[ 



and 



T > sech 2 



R < tanlr 



-dx 



I A' I 



-dx 



(65) 



(66) 



(67) 



The relevant transfer matrix is now 
E(x f ,Xi) = 








— 2iip 





dx 



(68) 



The Bogolubov coefficients are now 
E(oo, — oo) = 



a 13* 
(3 a* 



V exp 



-oo f' 







— 2iip 





dx 



(69) 



This type of analysis collects together and unifies several 
analytically known results for scattering over the barrier, 
such as for asymmetric square wells and Poschl- Teller po- 
tentials. (See Section VI.) After a few general comments, 
I shall turn to specializing this still rather general result 
to more specific cases. 



A. Reflection above the barrier 

The system (|||)-(|6^) can also be used as the basis of 
an approximation scheme for (3. Suppose that for all x 
we have \b{x)\ <C 1, so that \a(x)\ sa 1. Then 



5 



_ « _ exp(+^). 
This may be immediately integrated to yield 



-HOC:- 



exp (+2i<p) dx. 



(70) 



(71) 



J = Tin <J ^*^| = 2 Im{a&* exp(+2«/> tp )} • (78) 

For a properly normalized flux in the allowed region 
( | a | 2 — |6| 2 = 1), we have in the forbidden region 



2 Im {ab* exp(+2i<^ tp )} = 1. 



(79) 



Or the equivalent 



k'{x) 
k(x) 



exp 



-2i 



k(y)dy) dx. (72) 



This result serves to clarify the otherwise quite myste- 
rious discussion of "reflection above the barrier" given 
by Migdal |2^,|23| . Even though the WKB wavefunctions 
are buried in the representation of the wavefunction un- 
derlying the analysis leading to this approximation, the 
validity of this result for \(3\ does not require validity of 
the WKB approximation. 

If the shifted potential, V — Voo, is "small" we can 
recover the Born approximation in the usual manner. 



While this does imply 2|a||6| > 1, the inequality is unfor- 
tunately in the wrong direction to be useful for placing 
bounds on the transmission coefficient. 



C. Special Case 2-a 

Suppose now that V(x) is continuous and monotonic 
increasing or decreasing, varying from V-oo = V{— oo) to 
V+oa = V{+<x>). Suppose E > max-fV-oo, V +OQ } so there 
is no classical turning point. Then 



In that case k' 



exp(2« J k) i 
then yields 



mV'/{jrk) w mV'/(hrkoo), while 



-dx 



In 



^ + 00 



(80) 



exp(2ik oa x). A single integration by parts an d the transmission and reflection probabilities satisfy 



nrk 



(V{x)-V ao ) exp(+2ikoox) dx. (73) 



oo «/ — oo 



and 



B. Under the barrier? 

What goes wrong when we try to extend this analysis 
into the classically forbidden region? Analytically con- 
tinuing the system (p3|)-(|64|) is trivial, replace 



T > 



R < 



4k +no k- 



(fc+oc + fc— oo) 2 



(*4 



>) 2 



(&+oo ~\~ k— oo) 2 



(81) 



(82) 



(p'(x) = k^in = iy/2m(y - E)/h, 
and write 



ip(x) = iptp + i I n{y)dy, 
Jtp 



to obtain 



da k' 

— = +— b exp(-2z^ tp ) exp 

CLX £fx 

— = +— a exp(+2z< y 9 tp ) exp 

CLX £tx 



(74) 
(75) 

(76) 
(77) 



Thus we are violating our previous condition that ip be 
real, though we still require tp' ^ Q. This is a perfectly 
good Shabat-Zakharov system that works in the forbid- 
den region. But you cannot now use this to derive bounds 
on the transmission coefficient. The fly in the ointment 
resides in the fact that the formula for the probability 
current is modified, and that in the forbidden region the 
probability current is 



These bounds are immediately recognizable as the ex- 
act analytic results for a step-function potential Jl|,0,D, 
and the result asserts that for arbitrary smooth mono- 
tonic potentials the step function provides upper and 
lower bounds on the exact result. If we are interested 
in physical situations such as a time-dependent refrac- 
tive index |po] , p6| , or particle production due to the ex- 
pansion of the universe p7[ , this technique shows that 
sudden changes in refractive index or size of the universe 
provide a strict upper bound on particle production. 



D. Special Case 2-b 

Suppose now that V{x) has a single unique extremum 
(either a peak or a valley), and provided that E > 
max{I4o, ^cxtromum, V+oo} so that there is no classical 
turning point, then k(x) moves monotonically from fc_oo 



to k 



extremum 



and then back to fc. 



+oo • 



Under these circum- 



stances 



+ OC 



fc'l 



-dx 



"•extremum 



k- c 



in 



"--extremum 



extremum 



oc^+oc 



(83) 



(84) 



G 



This implies 

\a\ < cosh 
Which yields 

1/31 < sinh 



hi 



•'cxtremum 



In 



+00 



•-cxtrcmum 



To be more specihe, if in addition V(— 00) 
so that /c_oo = fc+oo, then we have 



and 



a < 



< 



k 2 4- P 

cxtrcmum ^ oo 



2k, 



cxtrcmum"' 00 



cxtrcmum 



(85) 
(86) 

= V(+oo), 
(87) 



2k, 



cxtrcmum^oo 



Translated into statements about the transmission and 
reflection probabilities this becomes 



T > 



and 



(E — Voo)(E — Vextremum) 



rcmum 



cxtrcmum 



(89) 



R < 



4 (Kxtrcmum 



[E Voo){E Vextremum) 4" 4 (Vextremum Vqo) 2 



Equivalently 



and 



T > 1 - 



R < 



(V e3 



y2iE Vextremum ^oo ) ^ 



(K 



cxtrcmum 



(2E - V c 



cxtrcmum 



(90) 



(91) 



(92) 



For low energies, these results are weaker than the 
bounds derived under Special Case 1, [(|46|), ([!?])] and 
[(P#|), J49j)], but have the advantage of requiring more 
selective information about the potential. For high ener- 



gies, 



h 2 (v c 



cxtrcmum 



2m 



Vnoldx 



(93) 



the present result (when it is applicable) leads to tighter 
bounds on the transmission and reflection coefficients. 

Numerous generalizations of these formulae are possi- 
ble. For example, at the cost of a little extra notation, we 
also already have enough information to provide a bound 



on an asymmetric barrier or asymmetric well, as long as 
it has only a single extremum (maximum or minimum) 
we apply the previous equations to derive 



and 



a < 



< 



k 2 

cxtrcmum 



1k c 



oo n — oo 



I*. 



2 

cxtrcmum 



^+oo ^— oo | 



(94) 



(95) 



Translated into statements about the transmission and 
reflection probabilities this becomes 



T > 



dk, k k 2 



and 



R< 



{^extremum 

Ik 2 

\ cxtrcmum 



&+00 00 } 
^+00 ^—00 } 



{^cxtrcmum ^+00^-00} 



(96) 



(97) 



Equivalently 



T > 



A{E — Vextremum) y/ (E ~ V +OQ ){E — V-oo) 
\{E ~ V extr emu m ) + V( E ~ V +00 )(E - Vloc)] 2 



(98) 



and 



R < 



[(E — Vextremum) ~ y/ (E — V +00 ){E — Vloo)] 2 
[{E - Vextremum) + y/iE^V+^E - VL^ft' 



(99) 



This can be compared, for example, with known analytic 
results for the asymmetric square well, see equation (144) 
in Section VI. 



E. Special Case 2-c 

Suppose now that V(x) has a number of ex- 
trema, (both peaks and valleys). I allow V(+oo) ^ 
V{— oo), but demand that for all extrema E > 
maxjVLoo, V+oo, Vex trcmum } so that there is no classical 
turning point. 

For definiteness, suppose the ordering is: — oo — > peak 
— > valley ... valley — > peak — > +00. Then 



+ OC 



k' 



dx = 



In 



In 



k 1 

^pcak 



1,1 

"'peak 
'Valley 



hi 



hi 



k 1 

"-valley 
peak 



v -(-oo 



peak 



(100) 



Defining 
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n p (fc) 

U v (k) 

n e (fc) 



nk i 
"-pcaki 

peaks 

\\_ ^valley ) 
valleys 

cxtrcn 



we see 



+°° I7„'l 



-dx 



In 



E|(A) 



fc— oak+oo^v (fc) 



This bounds the Bogolubov coefficients as 

iyfc) 



lal < cosh 



In 



n„(fc) 



That is 



and 



a < 



< 



fc_ 00 fc +00 ng(fc) + ng(fc) 
Ifc-oofc+oon^fc) - n2(fc)| 



(101) 
(102) 
(103) 

(104) 

(105) 
(106) 
(107) 



+oo ooiie (&) 

Then the transmission and reflection probabilities satisfy 

4/c +0 ofc-oon;?(fc) 



T > 



and 



R < 



{ng(fc) - k +00 k^ul{k)Y 
(n2(fc) + k+^k^uKk)} 2 



(108) 



(109) 



In these formulae, peaks and valleys can be interchanged 
in the obvious way, and by letting the initial or final peak 
sink down to V±oc as appropriate we obtain bounds for 
sequences such as: — oo — * valley — * peak ... valley — > 
peak — * +oo, or: — oo — > peak — > valley ... peak — * valley 
— > +oo. In the case of one or zero extrema these formulae 
reduce to the previously given results. [Equations (|96|)- 
(|97|).] Further modifications of these formulae are still 
possible, the cost is that more specific assumptions are 
needed to derive more specific results. 



V. PARAMETRIC OSCILLATIONS 

Though the discussion so far has been presented in 
terms of the spatial properties of the time-independent 
Schrodinger equation, the mathematical structure of 
parametrically excited oscillations is identical, needing 
only a few minor translations to be brought into the cur- 
rent form. For a parametrically excited oscillator we have 



dt 2 



u(t) 2 ^. 



(110) 



Just map t — > x, uj(t) — > k(x), and <fi — > tp. In the general 
analysis of equations (|28])-(|35|) the quantity i? should be 
replaced by 



% (t , UJ (t) = ± — . (Ill) 

2|<//| 

The analysis then parallels that of the Schrodinger equa- 
tion. Some key results are given below. 



A. Special Case 1 

If u)(— oo) = luq = uj(+oo) 0, then by choosing the 
auxiliary function to be tp = LUgt we can use equations 
[©-©] to deduce 



" : <osh(J— [ \u 2 (t) - u%\db I , (:i"l.2) 



\2uj 



and 



< sinh f-J- J \uj 2 {t) - uj 2 \dtj 



(113) 



B. Special Case 2 

If u>(— oo) and uj(+oo) ^ are both finite so that suit- 
able asymptotic states exist, and assuming ui 2 {t) > 
so that the frequency is always positive, then applying 
equations [(|66|)-(|67|)1 to the case of parametric resonance 
yields 



and 



< cosh 



< sinh 



+ OC 



K*)l 



dt 



dt 



(114) 



(115) 



C. Special case 2-a 

Suppose now that w 2 (i) is positive semidefinite, con- 
tinuous, and monotonic increasing or decreasing, varying 
from w_oo = uj{— oo) ^ to lu+oc = oj(+oo) ^ 0. The 
Bogolubov coefficients satisfy 



lal < 



CJ_oo + U> 



oo i ^+00 



and 



< 



2\/tJ-oo^+oo 



oo — ^+oo| 
2^/W-co^+oo 



(116) 



(117) 
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D. Special Case 2-b 

Under the restriction ui{— oo) = ujq = lu(+oo) 0, 
with the additional constraint that u(t) has a single 
unique extremum (either a maximum or a minimum but 
not both), and provided that <^ xtromum > so that we do 
not encounter complex frequencies (no classical turning 
point), the Bogolubov coefficients satisfy 



a < 



extremum 



extremum 



and 



< 



extremum I 



2 LdQ ^extremum 



(118) 



(119) 



Suppose now that uj 2 (t) has a single unique extremum 
(either a peak or a valley), but that lu(+oo) u>(—oc), 
and further that 0J 2 (t) > so that there is no classical 
turning point. The Bogolubov coefficients satisfy 



a < 



extremum 



2^/ U) — oo ^+ oo ^extremum 



and 



< 



^extremum 



(120) 



(121) 



E. Special case 2-c 

Suppose now that u>(t) has a number of extrema (both 
peaks and valleys). I allow lu(+oo) ^ lu(— oo), but de- 
mand that for all extrema Wg Xtremum > so that there is 
no classical turning point. 

For defmiteness, suppose the ordering is: — oo — ► peak 
— > valley ... valley — > peak — > +oo. Define 



HpM = Y[ W pcak, 
peaks 

n v (u) = J| < aUcy , 

valleys 



extremum* 



n e(^) = n 

extrema 

The Bogolubov coefficients satisfy 
and 



< 



^/W+ooW-ooITe^) 

(fc) -ng(fc) | 

^CJ+ooW-ooIL^fc) 



(122) 
(123) 
(124) 

(125) 
(126) 



In these formulae, peaks and valleys can be interchanged 
in the obvious way, and by letting the initial or final peak 
sink down to oj±oo as appropriate we obtain bounds for 
sequences such as: — oo — > valley — > peak ... valley — > 
peak — * +oo, or: — oo — * peak — > valley ... peak — * valley 
— > +oo. In the case of one or zero extrema these formulae 
reduce to the previously given results. 

Again, further specializations of these formulae are 
still possible. As always there is a trade-off between the 
strength of the result and its generality. 



VI. COMPARISON WITH KNOWN ANALYTIC 
RESULTS 

For comparison purposes, in this section I collect sev- 
eral known analytic results and show how they relate to 
the general results presented in this article. 



A. Delta-function potential 

For a delta function potential 

V(x)=a5(x), (127) 
the transmission coefficient is known to be |||| 

1 



T 



1 



ma" 
2Eh 



(128) 



This satisfies the bound (|4q), and also (|48|), and for E — > 
oo asymptotically approaches the bound, thus showing 
that the bound cannot be improved in the high-energy 
regime unless additional hypotheses are made. 

Though these bounds were derived assuming well- 
behaved functions the statements ( pro]) and (|^) continue 
to make good sense even for delta-function potentials. 
Thus any smooth set of well-behaved functions tending to 
a delta-function limit may be used to establish ( ff6| ) and 
( ptcf ) even for potentials containing delta- function contri- 
butions. 



B. Double-delta-function potential 

For the double delta function 

V{x) = a {5{x -L/2) + S(x + L/2)} , 
the transmission coefficient is |ll|] 

1 



T 



1 



Iraq . 



} cos(fcX) 



1 ( 2ma \ 

2 \~hFk) 



' sin(fcL) 



(129) 



(130) 



It is an easy exercise to check that this satisfies the 
bounds @ and ©. 
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C. Square barrier 

Tunneling over a square barrier is an elementary prob- 
lem which however is not always discussed in the text- 
books. (Tunneling under a square barrier is much more 
popular.) The exact transmission coefficient is 



E(E-V e 



(131) 



E(E - V e ) + \V 2 8m 2 (y/2m(E-V e )L/h) ' 
(See Landau-Lifshitz |l| or SchifT ||.) If we re- write this 



as 



T 



1 



1 



mV 2 L 2 sin 2 (y/2m(E-V e )L/h) 



(132) 



2Eh 2 



2m(E~V)L 2 /h 2 



then it is clear that the bound ( (46[ ) is satisfied. It is 
also possible to verify that this satisfies the general lower 
bound ( |66| ) that I have presented above, and in fact os- 
cillates between this lower bound and the upper T < 1 
unitarity limit. For certain values of the barrier width 
[fcextremum-k = (2n + l)7r/2] the square well saturates this 
bound thus showing that this bound cannot be improved 
unless additional hypotheses are made. 



proving this is tedious.) Start by noting that for this sech 
potential 



T > tanh 2 [7iV2?7j£X/?i], 
and use the inequality (x > 0) 



tanh x > 



l + x< 



> sech 2 (l/a;). 



Then 



T > sech 2 [h/(wV2mEL)] 



sech 



to 2LIVJ 



tt\2E h 8m\V e \L 2 



(137) 

(138) 

(139) 
(140) 



Provided that the extremum is a peak, V^cak > we can 
use the bound SmYvcakL 2 < ft 2 to deduce 



T > secli 



2i| Vpcakl 

2E h 



(141) 



This is the particularization of (|34|) to the present case. 
If V e < we need a different analysis. 



D. Tanh potential 

For a smoothed step function of the form 



V[x) = 



V+ °° 2 V ~°° tanh (I), (133) 



the reflection coefficient is known analytically to be 



R 



sinh[27r(fc_ 



sinh[27r(fc_oo + k +00 )L} 



(134) 



This certainly satisfies the general bounds (^l|)-(|8^) 
enunciated above, and as L — > approaches and satu- 
rates the bound. 



E. Sech potential 

For a sech 2 potential of the form 

V(x) = V e sech 2 (x/L), (135) 
the transmission coefficient is known analytically to be 

T _ sinh 2 [7rV '2mEL/K\ 

smh 2 [K^2mEL/h] + cos 2 [|7rv/l - 8mV e L 2 /h 2 ] 

(136) 

provided 8mV e L 2 < Ti 2 . This satisfies the general 
bounds, both (|3) and J8l|), enunciated above. (Though 



F. Asymmetric square-well potential 



For the asymmetric square well 




x < a; 

a < x < b; 

b < x. 



(142) 



we define ki 
ficient is |19 



y/2m(E — Vi)/h. The transmission coef- 



T : 



Akxk 2 ^ 

k 3 ) 2 k 2 + [k 2 k 2 + k 2 (k 2 - 



Then 



T > 



4fcifc|fc 3 



{k 2 + k!k 3 y 



k\ -ki)] sin 2 (k 2 L)' 
(143) 



(144) 



Similarly to the case for the symmetric square well, the 
transmission probability for the asymmetric square well 
oscillates between the bound (Q) and the unitarity limit 
T = 1. For certain values of the width of the well 
[k 2 L = (2n+l)7r/2] the transmission coefficient saturates 
the bound thus showing that this bound cannot be im- 
proved unless additional hypotheses are made. Because 
V-oo 7^ V+oo the bound (gj) is not applicable, at least 
not without modification from its original form. 
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G. Poschl- Teller potential 



ACKNOWLEDGMENTS 



For the Poschl- Teller potential 
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we have 

V-oo = V e- 2 ^, Vcxtrcmum = 0; V-oa = . 

(146) 

The transmission coefficient is 
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2 sinh(7rfc_ 00 L) sinh(7rfc+ 00 i) 



cosh[7r(fc_ 00 + k +oc )L] + cos 



(147) 



It is now a straightforward if tedious exercise to check 
this analytic result against all the bounds derived in this 
article. 



VII. DISCUSSION 

The various special cases discussed above are merely 
specific examples of the general results (|3^)-(|33|) and 
(|34|)-([35|) illustrating the power of the technique. There 
are many other variations on the bounds presented above 
that can be derived for specific choices of tp(x) and spe- 
cific restrictions on the scattering potential V(x). 

The most general form of the bounds are given in 
equations ©-© and (0)-©. Because of the large 
amount of freedom in choosing the function (p these 
bounds encode even more specific cases beyond those 
discussed in this article, and have the potential for lead- 
ing to new interesting specific cases. The special cases 
I discussed in this article were chosen for directness and 
simplicity. 

For instance, Special Case 1, as presented in equations 
(H)-(0) and has the advantage that it ap- 

plies to both scattering over the barrier and under the 
barrier. On the other hand, Special Case 2, as presented 
in equations (|66])-(|67|) and their specializations, applies 
only to scattering over the barrier but has the advantage 
of being much more selective in how much information is 
needed concerning the scattering potential. 

In summary, the bounds presented in this article are 
useful in establishing qualitative analytic properties of 
one-dimensional scattering, and as such are complemen- 
tary to both explicit numerical investigations and the 
guidance extracted from exact analytic solutions. 



[1 

[2. 
[3. 
[4 
[5, 

[<: 

[8 
[9 
[10 

[11 
[12 
[13 
[14 

[is; 

[16 

[17; 
[is; 

[19 

[20 



visser@kiwi.wustl.edu 

L.D. Landau and E.M. Lifshitz, Quantum Mechanics: 
Non-relativistic theory, (Pergamon, New York, 1977). 
G. Baym, Lectures on Quantum Mechanics, (Benjamin, 
New York, 1969). 

S. Gasiorowicz, Quantum Physics, (Wiley, New York, 
1996). 

E. Merzbacher, Quantum Mechanics, (Wiley, New York, 
1965). 

J. Singh, Quantum Mechanics: Fundamentals and appli- 
cations to technology, (Wiley, New York, 1997). 
P.M. Mathews an K. Venkatesan, A textbook of Quantum 
Mechanics, (McGraw-Hill, New York, 1978). 
A.Z. Capri, Non-relativistic Quantum Mechanics, 
(Benjamin-Cummings, Menlo Park, California, 1985). 
See esp. pp. 95-109. 

P. Stehle, Quantum Mechanics, (Holden-Day, San Fran- 
cisco, 1996). See esp pp. 57-60. 

L.I. Schiff, Quantum Mechanics, (McGraw-Hill, New 
York, 1955). 

C. Cohen-Tannoudji, B. Dui, and F. Laloe, Quantum Me- 
chanics, (Wiley, New York, 1977). 

A. Galindo and P. Pascual, Quantum Mechanics I, 
(Springer- Verlag, Berlin, 1990). 

D. Park, Introduction to the Quantum Theory, (McGraw- 
Hill, New York, 1974), 

A.T. Fromhold, Quantum mechanics for applied physics 
and engineering, (Academic, New York, 1981). 
M. Scharff, Elementary Quantum Mechanics, (Wiley, 
London, 1969). 

A. Messiah, Quantum Mechanics, (North-Holland, Ams- 
terdam, 1958). 

K. Chadan and P.C. Sabatier, Inverse problems in 
quantum scattering theory, (Springer- Verlag, New York, 
1989). 

W. Eckhaus and A. van Harten, The Inverse Scatter- 
ing Transformation and the theory of Solitons, (North- 
Holland, Amsterdam, 1981). 

R. Peirls, Surprises in Theoretical Physics, (Princeton, 
Princeton, 1979). See esp. pp. 21-22. 
M. Bordag, J. Lindig, and V.M. Mostepaneko, Class. 
Quantum. Grav. 15 581 (1998). 

N. Froman and P.O. Froman, JWKB approximation: 



11 



Contributions to the theory, (North-Holland, Amster- 
dam, 1965). 

[21] Transfer matrix techniques are discussed, at varying lev- 
els of detail, by Merzbacher §, Singh g, and Math- 
ews jf|. 

[22] A.B. Migdal, Qualitative methods in quantum theory, 

(Benjamin, London, 1977). 
[23] A.B. Migdal and V.P. Krainov, Approximation methods 

in quantum mechanics, (Benjamin, London, 1969). 
[24] P.M. Morse and H. Feshbach, Methods of Theoretical 

Physics, (McGraw-Hill, New York, 1953). 
[25] S. Liberati, M. Visser, F. Belgiorno, and D.W. 

Sciama, Sonoluminescence: Boqolubo v coefficients for the 

QED vacuum of a collapsing bubble, quant-ph / 980502c ; 

Sonolumine scence as a QED vacuum effect, [quant- 

ph/980503ll . 

[26] E. Yablonovitch, Phys. Rev. Lett. 62, 1742 (1989) 
[27] N.D. Birrell and P.C.W. Davies, Quantum fields in curved 
space, (Cambridge University Press, Cambridge, Eng- 
land, 1982). 



12 



